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ABSTRACT 

The reahzation that direct imaging of extrasolar planets could be technologically feasible 
within the next decade or so has inspired a great deal of recent research into high-contrast 
imaging. We ourselves have contributed several design ideas, all of which can be described as 
shaped pupil coronagraphs. In this paper, we offer a complete and unified survey of asymmetric 
shaped pupils designs, some of which have been published in our previous papers. We also 
introduce a promising new design, which we call barcode masks. These masks achieve the required 
contrast with a fairly large discovery zone and throughput but most importantly they are perhaps 
the easiest to manufacture and might therefore stand up best to a refined analysis based on vector 
propogation techniques. 



1. Introduction 

The discovery of more than 100 extrasolar Jupiter-sized planets in just the last decade has generated 
enormous interest, both among astronomers and the public, in the problem of discovering and characterizing 
Earthlike planets. NASA is already planning its next large space-based observatory, the Terrestrial Planet 
Finder (TPF)^ with a planned launch date toward the end of the next decade. TPF's primary objective will 
be to discover Earthlike planets and characterize them for indications of life. 

While the technical challenges for TPF are great, foremost among them is the problem of high-contrast 
imaging. In order to discover as many planets as possible, it is necessary to design an imaging system that 
achieves very high contrast between the parent star and the nearby planet. An earlier study by Brown et al. 
(2002) indicates that a D = 4m class visible-light (i.e. 400nm < A < 650nm) instrument ought to be able to 
discover about 50 extrasolar Earth-like planets if it can provide contrast of 10~^^ at an angular separation 
of 3X/D and that a 4 x 10m class telescope ought to be able to discover about 150 such planets if it can 
provide the same contrast at a separation of 4 A/I). 

Many approaches have been examined for designing a telescope with this level of contrast. The most 
promising fall into two broad categories — nulling interferometers (operating in the infrared) and coronagraphs 
(operating in the visible). One important subset of coronagraphs, referred to as shaped pupils, has been 
gaining interest (Spergel (2000); Kasdin et al. (2003); Vanderbei et al. (2003a,b)). These are apodized 
entrance pupils that rely solely on one/zero binary openings. In this paper we present a unified treatment of 
the one-dimensional shaped pupil designs, a few of which have been previously published, and we introduce 
a promising new design which we call barcode masks. 

The paper is organized as follows. In the next section we briefiy review the relationship between 
pupil-plane apodization and the corresponding image-plane electric field and point-spread function. In the 
following section we discuss the various performance metrics that characterize a planet finding telescope. 



- 2 - 



In Section 4 we show how optimization problems can be formulated and used to find the "best" entrance 
pupil apodization for high contrast in one- dimensional square apertures as well as azimuthally symmetric 
(circular or elliptical) pupils. In Section 5 we describe how these results can be used to find various families 
of assymmetric shaped pupils. In section 6 we present a preliminary sensitivity analysis to determine the 
tolerances of the pupil approaches to manufacturing errors. 

2. Apodization 

The main objective of our work is to develop optimal, high- contrast, binary shaped pupils (that is, pupils 
that consist of only open and closed areas). Since pupil masks can be viewed mathematically as a special 
case of apodized pupils (where the apodization function is simply zero-one valued), we begin by developing 
some notation and concepts in the general context of pupil apodization. In this paper we use scalar far-field 
diffraction theory (i.e., Fraunhoffer optics). We are currently investigating near field and vector propagation 
effects and will report on those more refined models in a future paper. 

In the Fraunhoffer regime, the image-plane electric field E{^X)^ produced by an on-axis point source 
and an apodization function^ A{x,y)^ reduces to the Fourier Transform of the apodization: 

^^^'^^ " / (1) 
where S denotes the overall aperture and can be taken, for example, to be a square: 

S = {{x, y) : -1/2 <x< 1/2, -1/2 <y< 1/2} 

Here, x and y denote the coordinates in the pupil plane measured in units of the aperture D and (^, 
Q denote angular (radian) deviations from on-axis measured in units of wavelength over aperture (A/D) or, 
equivalently, physical distance in the image plane measured in units of focal-length times wavelength over 
aperture (fX/D). 

When S is a square and the apodization function, A{x), is one-dimensional and symmetric about the 
i/-axis, the expression for the electric field can be simplified into a single, one-dimensional Fourier Transform: 

C) = ^ / A{x) cos{2^ix)dx (2) 

In this paper we will only be considering apodizations (and shaped pupils) along a single dimension. 
For both square and circular apertures, the point spread function (PSF) is the square of the electric field 
and, for planet finding, the apodization is chosen so that the corresponding PSF has the desired contrast. 

3. Performance Metrics 

The goal of our apodizations is to create point spread functions that provide the contrast necessary 
to detect a planet with as short an integration time as possible. In order to compare different apodization 
approaches and to perform optimizations it is therefore necessary to introduce a few important performance 
metrics (Vanderbei et al. (2003b)). The first is contrast^ which is a function of the linear position in the 
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image plane. In rectangular coordinates we define contrast by: 



e2(c,c)/£;'(o,o). 



We say that we have high- contrast when this ratio is very small. For planet finding, the consensus is that 
a contrast of 10~^^ outside a predefined radius is necessary in order to image an Earth-like planet that is 1 
AU from its Sun- like star. 

The second and third metrics are the inner working angle, ^iwa, and, of lesser importance, the outer 
working angle, ^owa- We wish to design pupil masks that have high-contrast for all ^'s in the interval 

The final performance metric is integration time, or its surrogate, throughput. Integration time is a 
critical metric for any planet finding telescope such as TPF. Smaller integration time relaxes requirements 
on system stability, simplifies control, allows more observations, and potentially allows return visits to 
interesting systems. In Kasdin et al. (2003) we described a number of different approaches to image analysis 
that resulted in quantifiable integration times that we used to compare different coronagraph approaches. 
All are based on some type of photometry on a potential planet and a predefined signal-to-noise ratio. 
While useful for comparing different coronagraphs, the resulting equations were too complicated to use in 
optimization work. Therefore, we consider here slightly simpler measures of throughput. (These measures 
are closely related to integration time because they measure the amount of light entering the system, but 
they do not account for the effect of PSF shape. While we acknowledge that PSF shape or "sharpness" is an 
important quantity in determining absolute integration times, it does not have a strong effect on the relative 
comparisons or optimizations.) Even asking only for measures of throughput is not such a simple question 
as there are multiple ways one might formulate an answer. 

Perhaps the most natural measure of throughput is the amount of energy that falls into the main lobe 
of the PSF relative to the total energy conveyed through a fully open aperture: 



where Ciwa ^iwa denote the locations in (^, Q of the first null of the PSF. We call this the Airy-throughput 
(Vanderbei et al. (2003b)). 

There are two other relevant measures of throughput. The total throughput of an apodization is the 
integral of the PSF over the entire image plane, which, by Parseval's theorem, is the integral of the square 
of the apodization over the entire entrance pupil. 

Finally, the electric field at (^ = 0, C = 0), 



for one-dimensional rectangular apodizations provides another measure of the "central" throughput of the 
apodization, since its square is the peak throughput density at the center of the Airy disk. If A{) is zero-one 
valued, then the apodization can be realized as a shaped pupil mask. In this case, ^(0, 0) is precisely the 
open area of the mask (and is also the total-throughput). For this reason, we call -£^(0, 0) the pseudo-area of 
the apodization, even when the apodization is not zero-one valued. 




(3) 



.1/2 

E(0,0) = 2 / A{x)d: 



(4) 
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Fig. 1. — Left The PSF for the one-dimensional prolate spheroidal apodized square aperture (Slepian 1965) 
of unit area plotted on a logarithmic scale with black areas 10~^^ below brightest. Right On-axis cross 
section of the PSF showing an inner working distance of AX/ a. The total throughput is 25%. 



4. Optimal Apodizations for High-Contrast 

In this section we discuss the general problem of designing apodized pupils for high contrast. This 
problem is not new; an excellent survery can be found in Jacquinot and Roizen-Dossier (1964) with more 
recent ideas in Indebetouw (1990) and Watson et al. (1991). Nisenson and Papaliolios (2001) describe their 
concept of using an apodized square aperture for high-contrast, where discovery is made along the diagonals. 
We describe below alternative optimal apodizations that lead naturally to various shaped pupil concepts. 



4.1. One-Dimensional Apodizations 

Based on Eq. (2) we ask for the optimal one-dimensional apodization for planet finding. Of course, there 
are a number of ways one can state the optimization problem, all involving shaping the point spread function 
(PSF) in some way. For instance, Jacquinot and Roizen-Dossier (1964) describe many different approaches 
to designing optimal apodizations. The most important optimal apodization is due to Slepian (1965), who 
introduced the prolate spheroidal wave function as the optimal 1-D telescope apodization. Slepian's one 
dimensional optimization problem asks for the function that concentrates as much light as possible into the 
central lobe of the finite Fourier transform. Slepian and Pollack (1961) derive a finite fourier transform 
analog of the uncertainty principal and show that the function that solves this optimization problem is the 
zero order prolate spheroidal wavefunction, that is, the solution to the wave equation in prolate spheroidal 
coordinates (there is also an elegant solution to this problem using the calculus of variations). Their statement 
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Fig. 2. — Left A one-dimensional barcode mask optimized for a dark region in [0,40] (all slots are not visible 
because of resolution limitations). Center The 2-D point spread function of the barcode mask Right A cross 
section of the optimal one-dimensional PSF near the inner working angle. The Airy throughput is 25 %. 



is equivalent to the computationally easier problem: 

POO 

minimize / E{^^0)'^d^ 
subject to ^(0) = 1. 

Choosing ^i^a = 4 provides the needed 10~^^ contrast for ^ > ^iwa- Fig. 1 shows the PSF resulting from a 
square aperture apodized in one dimension with a prolate spheroidal wave function. The total throughput 
of this apodized system is 25%. 

We note that outside the iwa the PSF has higher than needed contrast, which leads us to a slightly 
different optimization problem — namely, to maximize the pseudo-area subject to contrast constraints. This 
has the enormous advantage that it can be formulated as an infinite dimensional linear programming problem: 



maximize Jq A{x)dx 

subject to -10-^^(0, C) < ^(0, C) < 10-^^(0, C), 6wa < ^ < ^owa, (5) 

< A{x) < 1, 0<x < 1/2, 

Discretizing the sets of x's and ^'s and replacing the integrals with their Riemann sums, problem (5) is 
approximated by a finite dimensional linear programming problem, which can be solved to a high level of 
precision (see, e.g., Vanderbei (2001)). The numerical solution to this problem reveals that the optimal 
solution is zero-one valued. In other words, the optimal solution is, in fact, a binary mask consisting of a 
series of slots in the ^-direction of varying width alternating along the x-axis. We call such a mask a barcode 
mask. Fig. 2 shows a cross section of the PSF of such a mask optimized for high contrast in the region 
4 < ^ < 40 {\/D). The Airy throughput of this mask is 25%. Note that this response is very similar to 
the sine-squared bandlimited coronagraph, with bands of discovery zones in the image plane (Kuchner and 
Traub (2002)). 

In order to investigate as many planets as possible (particularly at larger distances), there is a strong 
desire to achieve as small an inner working angle as possible. One benefit of small inner working angle is a 
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reduction in the size of the telescope required to image the fuh habitable zone of these stars. One way to 
achieve smaller inner working angle is to optimize a narrower discovery region (that is, bring in the outer 
working angle). Fig. 3 shows a cross section of the PSF of such a design optimized for high contrast in the 
range 3 to 5 A/I) with an Airy throughput of 33.8%. It also possible to design a mask for an inner working 
angle as small as 2 A/D, but the width becomes prohibitively small (2 to 2.2 X/D). 

Perhaps a more natural, but more complicated, optimization is to maximize the total throughput rather 
than the pseudo-area. This results in a nonlinear optimization that is much more difficult to solve. Never- 
theless, we examined some solutions and found the same result as the pseudo-area problem. We also note the 
importance of considering the Airy throughput for these finite band masks rather than the total throughput 
as the light gathered from a planet with this PSF would correspond only to the inner core. It is thus more 
appropriate to ignore the light scattered outside the outer working distance. 

Finally, it is possible to numerically obtain a smooth apodization similar to the one in Fig. 1 by 
modifying the original optimization problem. To accomplish this, we add seemingly artificial smoothness 
constraints. Motivated by the fact that optimal apodizations look qualitatively like a Gaussian function, we 
impose smoothness constraints that correspond to the following conditions: 

log(A)' < (6) 
log(A)" < (7) 

The resulting max-pseudo-area optimization problem is: 

maximize J^^'^ A{x)dx 

subject to -10"^£;(0, C) < ^(0, C) < 10-^E{0, C), 

< A{x) < 1, 
A\x) < 0, 
A{x)A'^{x)<A'{x)^, 



iwa _ s — Sowa 9 

0<x< 1/2, (8) 
0<x< 1/2, 
0<x< 1/2. 



Fig. 4 shows the resulting smooth apodization function and the corresponding PSF (i.e., the Fourier Trans- 
form). Not surprisingly, it has a strong similarity to the prolate spheroidal wavefunction. 
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Fig. 4. — Left A smooth one-dimensional optimal apodized function. Right The corresponding cross section 
of the PSF showing an inner working distance of 4A/a and contrast better than 10~^^. 

4.2. Phase Shifts in Apodized Masks 

The central drawback of apodized pupils, and one of the significant factors that led us to consider binary 
masks, is the difficulty manufacturing them. It is extremely hard to create a smoothly varying mask to the 
necessary accuracy for 10~^^ contrast. In addition, there is no known process for creating apodized masks 
that doesn't also introduce a phase shift as a function of the transmission. This phase shift can dramatically 
reduce the performance of the apodized pupil. To better understand the issue and in an effort to ameliorate 
the problem, we examined optimizing in the presence of such a phase shift. 

Incorporating the effect of phase shift, the correct expression for the electric field is 

.1/2 
J-l/2 

where A: is a constant giving the rate that lag appears as attenuation is increased. Note that, unlike our 
earlier expressions, this complex expression does not reduce to a real- valued one. In Fig. 5, the upper 
apodization was computed with k = (i.e., no phase shift) and using ^iwa = 4 and ^owa = 40. The lower 
apodization corresponds to a nonzero phase shift. It was computed using a larger inner working angle of 
6wa = 4.5. The problem seems to be infeasible with a tighter inner working angle. 

If the apodization is computed assuming no phase shift but in reality there is a phase shift, then one 
expects the PSF to degrade because of the shift. Fig. 6 shows the amount of degradation one will see with 
k = 0.1. 

Finally, Fig. 7 shows the psf corresponding to the apodization computed assuming a phase shift of 
k = 0.1. Note that the iwa is now 4.5. 
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Fig. 5. — Optimal Apodization with and without Phase Shift 

5. Asymmetric Shaped Pupils 

Because of the problems associated with manufacturing and the limitations due to phase shifts described 
above, apodized pupils present unique and difficult challenges as a mechanism for modifying the point spread 
function to achieve high contrast. We favor, instead, shaped pupils, where manufacturing involves only 
cutting openings in the material at accuracies easily achieved with current technology. We have already, in 
fact, presented one class of shaped pupils as the optimal solution to the high-contrast imaging problem — the 
barcode mask of Section 4.1 (see also Vanderbei et al. (2003b) and Vanderbei et al. (2003a) for examples 
of circularly symmetric shaped pupil masks). In the following subsections we present a collection of other 
asymmetric shaped pupil masks derivable from a variety of optimization problems. 




5.1. Single Pupil 

David Spergel was the first to realize that shaped pupils could be used to achieve essentially arbitrarily 
deep nulls close to the central star. His first suggestion (Spergel (2000)) was to use a gaussian-shaped 
mask. Shortly thereafter, Jeremy Kasdin realized that Slepian's prolate-spheroidal wave function (Slepian 
and Pollack (1961)) could be used here too as the optimal shape for achieving high contrast along one axis 
of the image plane. This is seen by returning to the expressions for apodizations in a square or rectangular 
aperture (Eqs. 1 and 2). If in place of a smooth apodization we assume instead a single pupil opening of 
width w{x) relative to the ^ = axis we have for the electric field on the C = ^^is in the image plane: 

.1/2 

^(^,0)= / w{x)e-''^''^''dx (9) 

J-l/2 

In other words, the field on the axis of the image plane is just the one-dimensional Fourier Transform 
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Fig. 6. — PSF for Optimal Apodization without Phase Shift. Here is shown the PSF for the apodization 
computed without accounting for phase shift. The upper curve is the PSF one gets with the phase shift term 
using A: = 0.1. 

of the pupil width. Thus, an optimization problem for that width is identical to the Slepian problem we 
described earlier for the one-dimensional apodization. A single pupil with the Slepian shape would provide 
the needed contrast. Such a pupil is shown in Fig. 8. It achieves a contrast of 10"-^^ everywhere along the 
X-axis except within an Airy disk of radius 4A/D. This was our first proposed shaped pupil for TPF. 

It is clear from Fig. 8 that the disadvantage of the Spergel-Kasdin mask is the extreme narrowness of 
the high-contrast region at close working angles from the central star. Many rotations of the telescope would 
be necessary to image the entire discovery space, dramatically impacting integration time and mission life. 
The solution to opening up this discovery space is to introduce multiple openings, what we call multi-pupil 
masks. In fact, the barcode mask introduced earlier is an example of a multi-pupil mask. In what follows 
we describe several other approaches to accomplishing this end. 

5.2. Multi-Pupil Square Aperture 

The simplest possible multi-pupil mask involves just repeating the single pupil, one on top of another, in 
a square aperture until the desired PSF is achieved. In this approach, rather than optimizing a design directly, 
we attempt to reproduce, as closely as possible, the PSF of the one-dimensional apodization described in 
Section 4.1 and shown in Fig. 1. To begin, we assume that each opening in the mask is identical and stacked 
vertically across the pupil. The opening, y{x), in the upper quadrant of a 27V pupil mask, can then be 
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Fig. 7. — PSF for Optimal Apodization with Phase Shift. Here is shown the PSF for the apodization 
computed with accounting for phase shift. Inner working angle is 4.5A/£>. 



described by the set S given by: 

S = 
Yu = 

while for the lower pupils, y{x) G Yi, 



{{x,y) ■.-l/2<x<l/2,y{x)eY^} 

N-l 

u 



n=0 



D n w{x) D n 
~2N ^ 2 '2" AT 



w{x) 



N-l 



n=0 



D n + 1 w{x) D n w{x) 



(10) 

(11) 

(12) 



where w{x) denotes the half width of the pupil and D is the horizontal width of the aperture. For what 
follows, all lengths are normalized by D to simplify notation. Note that the shape of each pupil is the 
same and determined by w{x). Again, the goal is to determine choices for this function that will yield an 
image-plane PSF matching the PSF corresponding to a given apodization (such as that in Fig. 1). The 
electric field for such a multi-opening mask can by found from Eq. (2): 



2s in(7rC) r'^'y^' 



1/2 N-l 



sin 27r( 



n + 1 w{x) 
~2N 2~ 



— sin 27r( 



n w(x) 

\ — 

2N 2 



cos(27r^x) 
sin(7rC) 



dx (13) 



which can be simplified by completing the finite sum: 

ri/2 



2sin«) r^/' / , 



Cw{x)) - ^^^"1^1^^^^^ sin«^(x)) ) cos{27r^x)dx (14) 



This looks like the field due to an arbitrary apodization A{x) in Eq. (2), except with A{x) replaced by 
the expression in parentheses. If we expand this expression in a Taylor series about = 0, we find that the 
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Fig. 8. — Top. The Spergel-Kasdin prolate-spheroidal mask. Bottom. The associated 2-D point spread 
function and its x-axis slice shown in decibels (10~^^ = — lOOdB). 

leading term in the series is 1 — 2Nw{x), independent of Thus, if we choose: 

^{x) = — (15) 

we can create an approximation to the continuous one-dimensional apodization response. What remains is 
to show that the higher order terms, for large enough TV, only produce interfering light outside some outer 
working angle Cowa- Replacing for w{x) from Eq. (15) into Eq. (14) results in: 

ttC Jo \ sm(7rC/TV) / 



We can now Taylor expand the term in parentheses in Eq. (16) for small {ir^/N) to find: 

E{^, = ^^^^ 1^^^ (^Mx) + O ' j cos(27rCa;)dx (17) 

Thus, for large enough N, the scattered light due to the finite number of openings can be made arbitrarily 
small within some range < C ^ Cowa- Fig. 9 shows the PSF resulting from such a multi-opening square 
aperture with 10 and 100 openings using the prolate spheroidal apodization shown in Fig. 1 for A{x). 

While it is certainly feasible, and perhaps not overly difficult, to manufacture a mask with over 100 
openings, it may be desirable to use fewer pupils. It is thus sensible to attempt to improve the resopnse with 
fewer openings by relaxing the requirement that each edge of the openings be the same, hoping that these 
extra degrees of freedom can be used to improve the match to the apodized pupil, albeit at the cost of a 
much more difficult computational problem. 

In this case, the set Yu for the upper half of the pupil is replaced by: 

N-l 

Yu= U 

n=0 



In bn{x) ln + 1 tn{x) 
2N ^ 2 '2^ 2~ 



(18) 



Fig. 9. — PSFs for a multipupil square aperture associated with the apodization shown in Fig. 1. Left A 10 
opening pupil. Right A 100 opening pupil. 



where bn{x) is the bottom edge of opening n and tn{x) is the top edge of opening n. A similar set holds for 
the lower half of the pupil. The resulting electric field in the image plane is given by: 



m,c) = 



2 sin(7rC) 



I 5 



siuTrC 



1 



N 



tn{x) 



sinK(^ 



cos(27r^x) 
sin(7rC) 



dx (19) 



Comparing again to Eq. (2) we find that making the correspondence: 



N- 



Aix) = J2 



1 SinTrC - tn{x)) - SinTrC (f + 6„(x)) 



n=0 



siuTrC 



(20) 



results in a multi-pupil mask reproducing the one-dimensional apodization. We again expand the right-hand 
side of Eq. (20) about C = and match the lowest order term (independent of Q: 



N-l 



N-1 



A{x) = l-^tn-^hr, 



n=0 



Unlike the previous multi-pupil mask, we can now use the extra degrees of freedom to make the coefficients 
of the terms of higher order in vanish. Each additional pupil pair allows us to eliminate two higher order 
terms. The computational problem to find the shapes of the openings becomes the solution to N simultaneous 
polynomial equations. While intuitively one would expect such extra degrees of freedom to result in a better 
match to the apodized solution with fewer pupils, the surprising result is that no significant improvement is 
made. Fig. 10 shows an example of a 6 opening pupil designed with this method, its PSF, and the PSF of 
a uniform 6 opening pupil. 

One of the advantages of asymmetric masks is that they can be made to have higher throughput by 
reducing the size of the high contast region (smaller opening angle). While this may be unattractive for planet 
detection, as many rotations become necessary to cover the discovery region, it can have large advantages 
for characterizing a planet at a known location. We summarize another category of assymetric masks ideally 
suited to this problem in the next section. 



Fig. 10. — Left An unequal six-opening pupil inscribed in a square aperture. Center The Corresponding PSF 
plotted on a logarithmic scale with black areas 10~^^ below brightest. Right The PSF of a 6 equal opening 
multi-pupil aperture. Spatial Frequencies shown from -20 to 20 X/D. 

5.3. Multi-Pupil Circular and Elliptical Apertures 

The multi-pupil asymmetric masks above are all embedded in square or rectangular apertures. Thus, 
the central core of the PSF is also square, resulting in slight degradation of the inner working angle off axis. 
Alternatively, multi-pupil assymetric masks can be designed that are embedded in circular and elliptical 
apertures by directly optimizing the edges of each opening relative to a performance function similar to 
those described in Section 4. Many such pupils are described in detail in Kasdin et al. (2003). Fig. 11 shows 
an example of a multi-openning pupil mask that was designed for an elliptical aperture. In an elliptical 
aperture, the asymmetry is used to leverage the long axis for improved iwa along the corresponding axis in 
the image plane. 




Fig. 11. — Left. A multiopening pupil mask designed to open up the high-contrast region. Right. The 
corresponding PSF. On the x-axis, contrast of 10~^^ extends over p > 4. 
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6. Tolerances 

A critical question for any apodized or shaped pupil coronagraph is how well the mask needs to be 
made. In this section we summarize the tolerancing requirements on an apodized pupil coronagraph and a 
barcode mask, insisting that any error produce scattered light below the 10~^^ contrast level. 

6.1. Apodized masks 

The electric field in the image plane for an apodized entrance pupil with some arbitrary error in ampli- 
tude and phase is given by a modification of Eq. (1): 

E{^X)= J J e^^''^y^A{x,y)e-^'''^''^^y^Uxdy (21) 

where 7(x, y) = y) + y) is an arbitrary complex function of position in the pupil representing the 
phase and amplitude transmission error of the mask. For simplicity, we consider here only the amplitude 
error in the apodization; we briefly discussed phase errors in Section 4.2. While the small amplitude error 
is arbitrary, because of the finite extent of the pupil we can model the error as a periodic function with 
fundamental period on [—1/2,1/2] (absorbing the top hat function into the apodization A{x,y)). This 
allows us to decompose the error into a Fourier Series: 

oo 

^a(x,y) ^^_^^J2an C0S(4 ' ^ + ^^n) (22) 
n=l 

where x = (x^y) is the position vector, kn = {kx,ky) is the wave vector of the harmonic term (with \kn\ 
being the spatial frequency). An and (pn are the amplitude and phase of each harmonic component, and e 
is small. This expansion allows us to write the electric field in the image plane as the nominal, designed 
for point spread function plus a small perturbation, E{^X) = ^o{^X) + ^^{^X)- Since AE is small, the 
intensity in the image plane is approximately given by: 

m, C)|2 « \Eo{^, C)|2 + 2\Eo{^, 0\\AE{^, 01 (23) 

Since the electric field in Eq. (23) is linear in the small change in intensity, superposition holds and it 
suffices to consider only a single harmonic component of the apodization error in Eq. (22). Thus, the field 
error is given by: 

^E = eJ J AnCOs{kn • X ^ ipn)A{x,y)e-'^'''^''^^y^^dxdy (24) 

By the Fourier Shift theorem, this is just a copy of the point spread function to the location {kx, ky) in the 
image plane, commonly referred to as speckle: 

AE = eEQ{£,-kx,C-K) 

The intensity error in the high contrast region, due to errors in manufacturing the apodization, is thus given 
by: 

A^^2e|i;o(e,C)l|i^o(C-fc.,C-fcj,)| <10-^°|i;o(0,0)|2 (25) 

The requirement that be less than 10~^° of the nominal PSF (below the contrast level) allows us 
to set a requirement on e, the magnitude of the apodization variation. Since we are concerned only with 



-15- 



the error in the high contrast region of the image plane, we can set \Eq{^, ()\/\Eo{0, 0)| to 10 ^ in Eq. (25), 
leaving: 

5x10-6 1^0 (0,0) I 

At the worst case error location of the speckle center, the ratio in Eq. (26) is unity, leaving a requirement 
on the manufacturing error of e < 5 x 10~^. Achieving a throughput variation with this accuracy is a 
formidable challenge. This, in part, motivates our search for shaped pupil solutions. 



6.2. Barcode Masks 



While it is demonstrably easier to manufacture a shaped pupil mask than an apodized one, it is still 
necessary to estimate the precision required. In this section we calculate the performance sensitivity of the 
barcode mask to manufacturing errors. Here we assume that the edge locations of each slot along the x-axis, 
r^, are in error by some function eifi{y) on the right half and e[f[{y) on the left half of the mask, where 
again e is small. We further assume that each ei and e[ are independent and identically distributed (i.i.d.) 
random variables with (^{e^} = and S'{eiej} = a'^Sij. 

The error in the electric field is then given by: 

iV 1/2 

= E / [(^^/^(^) - <fi^y^^ sin(2^r,0 + i{e,My) + e[fl{y)) cos(2^r,0] dy (27) 

i=l -^-1/2 

Here, rather than use Eq. (23) for the image plane error we instead turn to a stochastic analysis. Since 
(^{AE} is zero, we compute the expected value of the intensity: 

^{m, 01'} = \Eo{^, 01' + ^{\AE{^, C)|2} 



From Eq. (27) we therefore find. 



N 



<^{\AE{i,C)\'} = '^'Yl 



i=l 



1/2 



1/2 



+ 



J -1/2 



(28) 



< 10--'^|^o(0,0)|^ (29) 



As we did with the apodization error, we can model each of the mask edge error functions, fi{y)^ via a 
harmonic expansion. Since the small amplitude of the error was absorbed into the e^, the terms in brackets 
are bounded by one at the worst locations in the image plane. Thus, a conservative bound on the average 
edge error is given by: 

10-^1^0(0,0)1 



cr < 



/27V 



(30) 



For the barcode mask shown in Section 4.1, which has 44 edges per side, this translates into a manufac- 
turing accuracy on each edge of the mask of approximately 5 x 10" For a 2 inch mask, this corresponds to a 
25 nm accuracy requirement on the mask edges, something achievable with current fabrication technologies. 



7. Final Remarks 



In this paper we presented all the asymmetric shaped pupil masks we have studied to date for achieving 
high contrast imaging. All of these masks use apodization at a single pupil to shape the point spread function 
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of the imaging system, achieving 10~ contrast outside some desired inner working angle. These masks can 
be used in either rectangular or circular apertures, depending upon the specific mission tradeoffs. While the 
idea of using pupil apodization has been around for many years, the problems associated with manufacturing 
such smoothly varying masks to the accuracy needed for planet finding has yet to be solved. We show here 
that any apodized pupil can be made into a shaped pupil mask at a small cost in discovery space but 
with enormous practical benefit. Shaped pupils are far easier and less expensive to manufacture, are more 
tolerant to errors, and, in some cases, provide more throughput. While classical Lyot coronagraphs typically 
have more throughput, they suffer from the same manufacturing difficulties and are far more sensitive to 
alignments and pointing error. 
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